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The equation of state for dense fluids has been derived within 
the framework of the Sutherland and Katz potential models. The 
equation quantitatively agrees with experimental data on the isother- 
mal compression of water under extrapolation into the high pres- 
sure region. It establishes an explicit relationship between the 
thermodynamic experimental data and the effective parameters of 
the molecular potential. 



1. Introduction 

Deriving the equation of state (EoS) for water in a wide 
range of pressures and temperatures remains a challeng- 
ing open problem, especially in the high pressure region. 
In the recent papers [1-4] devoted to the water EoS, this 
problem was solved by fitting multiparameter formulae 
with such a large number of adjustable parameters that 
it approaches the number of experimental points. These 
methods are not based on reliable statistical mechanic 
foundations, and the applicability of these EoS is re- 
stricted. If the functional form of the EoS and their 
parameters are applicable to other substance or solution 
is an open question. 

In comparison with the majority of one-component liq- 
uids, water reveals many unusual properties in its nor- 
mal and supercooled states. The analysis of the diffusion 
peak of the quasi-elastic incoherent neutron scattering 
and the kinematic shear viscosity of water has shown 
that the global H-bond network disintegrates into an en- 
semble of weakly interacting clusters: dimers, trimers, 
tetramers, etc. [5-9]. It was also shown [10] that the 
properties of water in the supercritical region are deter- 
mined by the averaged spherically symmetric potential. 
Therefore, it is reasonable to use such well-known mod- 
els as Lennard- Jones, Buckingham, Sutherland or Katz 
potentials for water in the high pressure region. 



In this paper, we derive the EoS for a supercritical 
fluid within the framework of the Sutherland and Katz 
potentials. In our previous work [11], we used a new ver- 
sion of the thermodynamic perturbation theory (TPT) 
originally proposed by Sysoev [11[. The main feature of 
the proposed TPT is in the assumption that the func- 
tional form of the perturbed potential is identical to the 
potential of the reference system. Therefore, the devi- 
ation of the potential of the more compressed system 
from the potential of the less compressed system is con- 
sidered as a perturbation. On this basis, the concept 
of a reference thermodynamic state has been developed. 
A functional expansion of the free energy gave the pos- 
sibility to derive, at a certain choice of the parameter 
expansion, two EoS modifications within the framework 
of a realistic model and the "soft" sphere potential one. 
As was shown in [11], these EoS correctly described the 
isothermal compression for supercritical fluids of inert 
gases. 

Following [11[, we use the free energy perturbation 
expansion 
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Equation (1) can be transformed into the expression 
for pressure (the details are given in [11]) 
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Equation (6) can be expressed in the terms of 11 = 
y , using of the approximate quahty 
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11 = A+{1 + A) (L + K), 
where 

L = -^J g2ir)f{r)v{r)df, 



(7) 
(8) 



This expression was obtained within the framework of 
a reahstic potential model, that can be presented in the 
general form 



K=-^J g2{r)Av{r)df. 



(9) 



(^ (r) = $ (r) + -(/; (r) 



(3) 



where $ (r) is the repulsive part of the potential, and 
tp (r) is its attractive part. 

2. The Equation of State within the Framework 
of Sutherland and Katz Potential Models 

In the case of short-range potentials, the expression for 
the pressure can be rewritten in the form 

P = — ^ / .92 {r) f {r)v (r) dr- 

V 



Now the problem of deriving the EoS reduces to the 
evaluation of integrals (8) and (9). First, we consider 
the Sutherland potential 
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do is the molecular diameter, and the potential well 
depth £ is defined by £ = -pr ■ Then ip (qr) is written 
as 
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where j [r) = e — 1, v [r) = —\/ip(qr)r 

is the virial of intermolecular forces and vo (r) = 
—Vif (qr) f Iq^i = — V(/? (r) f is the virial of intermolec- 
ular forces in the reference state. With regard for the 
expression for the pressure in the reference state 
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we rewrite relation (4) as 
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where d — ^/^do. We calculate the singular force 



dip{r) 
dr 



following [13] and obtain 
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Av (f) = -kT [S{r-d)-S{r- do)] r+ 



m/3 
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6 {r — d) is the Dirac delta function and Q {r ~ d) is the 
Heaviside step function. The expression for L takes the 
form 
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(14) 



where / (r 4, d) = lini / (r) [13]. The integral in (14) 

can be evaluated if we assume that g2{'r) ~ erl^^^"^) _ 
Then, in view of ^ rs 1 + A, we have 
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Here, 7 (x, y) is the upper incomplete gamma function 
with X = ™,7 , v = —4m-- Using the same assumptions 
for integral (9), we obtain 
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where (72 (^ i c?o) = 1™ 92 (t")- Since the parameter A 

1 — >d+0' 

is small, we can expand 52 {d) at a point do in the Taylor 
series 
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r is the term which comprises the incomplete gamma 
function 
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For the Sutherland model, it is expressed in terms of 
the second virial coefficient B2 (T) 
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However, the analysis of the experimental data for 
some dense fluids (water, argon, neon, krypton) re- 
vealed that this term can be neglected under the con- 
dition H w 10^ with an accuracy of 1%. It is the 
fairly wide range of thermodynamic variables, where the 
isothermal compressibility is low (A ^ 1) correspond- 
ing to the pressure interval 100-2200 MPa. The terms 
2{1 + C (T)) A^-{l + C (T)) A can also be ignored with 
the same accuracy. Finally, we arrive at the EoS within 
the framework of the Sutherland model 



(21) 



P = Po + [B (T) + Pa] e^iryT _ 1 ^ (j) 



where D {T) = I + A{2 + C{T)). Expression (21) is 
a three-constant equation of state with the adjustable 
parameters B{T), D{T), and A{T) (see Table 1). The 
parameter B{T) depends on the temperature. It is re- 
lated to the pressure caused by attractive forces 

B (T) = ^kTg2 (r | do) dg - Po = 



g2{d)-g2{rido)[l + C(T)A 
where C(T) - ^i^iMM 



(17) 



the function of the 
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reference state. Since, in the approximation A ^ 1, 
/ (r I rf) = exp (^^) - 1 and ^ < 1, Eq. (7) in the 
case of small compressibility takes the form 



P = Po + [B(T) + Po] 



eMTW- 1 ) [1 + A (2 + C (P))] - 



Table 1. 


The values of the EoS parameters 




T, K 


300 


350 


400 


450 


500 


550 


B{T), MPa 


9.63 8.55 7.70 6.55 5.10 


3.35 


C{T) 


2.76 2.65 2.59 2.45 2.39 


2.10 


A(T), K-i 


0.0042 0.0037 0.0032 0.0027 0.0024 


0.002 



= P^ - Po = PS - 



NkT 



(22) 



where Pq is the pressure caused by the repulsive forces in 
the reference state, and Pq is a part of the pressure due 
to the attractive forces. The parameter A{T), within the 
framework of the Sutherland model, is expressed by the 
formula 



A{T) = ^^^. 
rac 

If we use the Katz potential model 



$(f) = 
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(24) 
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Fig. 1. Processing of the experimental PVT-data for supercritical 
water 

where Rq is the mteraction radius of the attractive 
forces, ^^ is the interaction constant, Eq.(21) retains 
its functional form, but the parameter A{T) is defined 
by the formula A{T) = ''"''"7jir°''°^'' - 

The parameters B{T) and C{T) are defined in a sim- 
ilar way. 

3. Experimental PVT-data Analysis 

We used the technique from the previous paper [11] to 
process the PVT data and evaluate the EoS parameters. 
It turns out that Eq. (21) yields good agreement with the 
experimental PVT data [14] under the extrapolation to 
the high pressure region. The results of the comparison 
are presented in Figure. 

Of special interest are the values of the parameter 
A{T) , because it is related to the parameters of the po- 
tential, formulae (23) and (25) for the Sutherland and 
Katz models. On the basis of (23), we evaluated the val- 
ues of the potential well depth e at a fixed value m — 12 
commonly used in the Lennard- Jones model (Table 2). 

If we fix the value of e = 0.650 kJ/mole, which cor- 
responds to the SPC/E model, we obtain a variation 
of the softness parameter m with the temperature. In 
general, fixing e at the values for well-known poten- 
tial models such as SPC/Fw [15], TIP3P/Fw ]15], and 



Table 2. 
m = 12 in the 


The values of e (Sutherland potential) at 
temperature interval 300—550 K 


T, K 


300 


350 


400 


450 


500 


550 


£, JlO^/mole 0.493 0.562 0.650 0.769 0.866 1.038 



TIP5P/Ew [16] leads to a variation of the softness pa- 
rameter m with the temperature within the framework 
of the Sutherland potential. 

4. Conclusion 

The approach developed on the basis of the free energy 
perturbation expansion and a new version of TPT re- 
sulted in some universality for the EOS statistical foun- 
dation of the low weight molecular supercritical fluids. 
In principle, the concept of the thermodynamic refer- 
ence state implies that an initial state {Pq, Vq) on the 
isotherm corresponds to the reference system with the 
unperturbed potential, and every subsequent point on 
the isotherm (Pi, V2), . . . , (P„, Vn) corresponds to the 
system with the perturbed potential at the isothermal 
compression of the system. This modification of TPT 
allowed obtaining the EoS which exhibits good results 
under the extrapolation to the high-pressure region and, 
most importantly, establishes a relationship between the 
parameters of the model potential and the thermody- 
namic properties of substances. This relationship gives 
estimations for the values of the parameter e (Table 2). 
Interestingly, the values are of the same order of mag- 
nitude as the values of many well-known water mod- 
els (SPC/E, SPC/Fw, TIP3P/FW, TIP5P/Ew). How- 
ever, the temperature dependence of e in this region of 
thermodynamic variables indicates that the form of the 
Sutherland model is unsuitable for this high pressure re- 
gion. Nevertheless, these data can be used as additional 
information for calibrating the potential parameters in 
simulations. 
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PIBH5IHH51 CTAHY BO/IH B OBJIACTI MAJTHX 
CTHCJIHBOCTEH 

B.K). BapdiK, fl.A. Hepyx, G.B. UaeAoe, LB. XmaMioK 

P e 3 lo M e 

0;i,ep»caHO CTaTHCTHMHO o6rpyHTOBaHe pIbhhhhh CTany rycTHx 
cj5jiiOi/;iB y paMKax MOflejiefl noTeHii,iajiy CiosepjieHfla Ta KaLi,a. 3a- 
nponOHOBane pIbhhhhh CTaHy 3 bhcokoio tomhIctio y3rofl>KyeTBC5i 
3 eKcnepHMeHTajiBHiiMH ^anHMH no i30TepMiMH0My cthchchhio bo- 
^,11 npH eKCTpanojiHii,ii b o6jiacTb bhcokhx thckIb. BcTaHOBjieno 

KijIbKiCHHH 3b'5I30K mI^C e4>eKTHBHHMH napaMeTpaMH MOflejIbHHX 

noTeHLi;iajiiB 3 napaMCTpaMH pIbhhhhh CTaHy. 
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